Signals of additional Z' boson in e + e — > W + W at the ILC 
with polarized beams 



B. Ananthanarayan 
Monalisa Patra 

Centre for High Energy Physics 
Indian Institute of Science 
Bangalore 560 012, India 

P. Poulose 

Department of Physics 
Indian Institute of Technology Guwahati 
Assam 781 039, India 



Abstract 

We consider the possibility of fingerprinting the presence of heavy 
additional Z' bosons that arise naturally in extensions of the standard 
model such as E§ models and left-right symmetric models, through their 
mixing with the standard model Z boson. By considering a class of ob- 
servables including total cross sections, energy distributions and angular 
distributions of decay leptons we find significant deviation from the stan- 
dard model predictions for these quantities with right-handed electrons 
and left-handed positrons at ^=800 GeV. The deviations being less pro- 
nounced at smaller centre of mass energies as the models are already 
tightly constrained. Our work suggests that the ILC should have a strong 
beam polarization physics program particularly with these configurations. 
On the other hand, a forward backward asymmetry and lepton fraction 
in the backward direction are more sensitive to new physics with realis- 
tic polarization due to interesting interplay with the neutrino t- channel 
diagram. This process complements the study of fermion pair produc- 
tion processes that have been considered for discrimination between these 
models. 



1 Introduction 



The International Linear Collider (LLC) is a proposed high energy, high luminos- 
ity electron-positron collider with the mission of studying the standard model 
(SM) at high precision and to look for signals beyond the standard model pQ. 
It has been proposed that an initial beam polarization program can also signif- 
icantly enhance its capabilities in meeting these objectives, see ref. [2J. 

One of the important processes that will be studied at high precision at ILC 
with and without beam polarization [2] is W^-pair production. Phenomenological 
studies of this process within the SM and some extensions have been carried 
out in great detail 3, 4 starting many years ago. Since properties of the weak 
gauge bosons are closely linked to electroweak symmetry breaking (EWSB) and 
the structure of the gauge sector in general, detailed study of W physics will 
throw light on what lies beyond the SM. 

On the other hand, it is entirely likely that there are additional Z bosons, 
denoted by Z' in the TeV range, see for instance the review section in ref. [5]. 
These are present in several economical extensions of the standard model. With 
this strong motivation, signatures of such a gauge boson is searched for in the 
past and existing colliders. Direct and indirect searches at LEP as well as at 
TeVatron and other existing facilities provide bounds on the masses of this 
particle and on other model parameters. Direct searches at TeVatron put lower 
limits of 630 - 1030 GeV [7] and LEP 2 put limits of 673 - 1787 GeV [| , while 
electroweak precision analysis of LEP provide lower limits of 475 — 1500 GeV [Hj 
on the mass, depending on the model considered. Even if not directly produced, 
they can be finger printed easily as they would mix with the traditional Z a of 
the electroweak model. Thus VK-pair production process has a winning edge 
compared to fermion pair production when it comes to the effect of this mixing. 
This is because, as these Z' do not interact with the standard W bosons, the 
W^-pair production process is insensitive to the presence of Z 1 in the absence of 
mixing. In contrast, fermion pair production process is sensitive to the presence 
of Z' even in the absence of mixing. 

We note that although such mixing is highly constrained by precision mea- 
surements at LEP and by other existing experimental data [5], with the high 
statistics expected at the ILC for W-pair production, even such small mixing 
can be probed effectively. The new effects could be manifested in departures 
from the standard model cross section for W-pair production, and in various 
differential distributions and asymmetries. 

Recently in the context of the littlest Higgs model(LHM), which also con- 
tain Z' bosons, we demonstrated the utility of several simple distributions in 
fingerprinting the model [TO]. As a result, it may be worthwhile considering how 
other popular Z 1 models arising in Eq unification and so called left-right sym- 
metric models (LRSM) and alternative left-right symmetric models (ALRSM), 
which have also been considered by the CLIC Physics Working Group [11] , can 
subject themselves to a diagnosis. This is the main aim of the present work. 
We have studied the process at reference energies of 500, 800 and 1000 GeV, 
and find that effects are pronounced only at the higher energy. In addition, we 
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conclude that initial beam polarization can significantly enhance the diagnostic 
ability of the ILC. In particular, for the class of models we are considering these 
effects are of importance for right-handed electrons and left-handed positrons 
in the process we are considering. This may be understood as arising from the 
dominant i-channel SM contribution which gets erased away with this choice of 
beam polarization revealing the new physics contributing through the s-channel. 
The above argument holds in the case for realistic degrees of polarization of the 
beams; it is expected that about 90% electron polarization would be achievable 
along with a positron polarization of 60% [2J. We present our results for both 
these cases and find that there are interesting effects even for the latter case 
as the t- channel contribution which now survives, can play an interesting and 
effective role. 

We will use the observables considered by us in our earlier work, ref. [TO] , 
which are total cross sections, single energy distribution of the secondary lep- 
ton, lepton angular distribution and forward-backward (FB) and left-right (LR) 
asymmetries. In addition to the above we consider an important and useful 
energy-energy correlation of the type first considered in the context of some 
anomalous gauge couplings by Dicus and Kallianpur [12]. Despite its obvious 
utility it has not received much attention, and we will demonstrate how this 
correlation in combination with beam polarizations can extract important in- 
formation on Z' models. 

The scheme of this paper is the following: In Sec. [5] we discuss the details 
of the models we are considering. In Sec. [3] we consider the kinematics of the 
process and the subsequent decays in great detail. The section is organized in 
several subsections for convenience. In Sec. |4]we present a discussion including 
a comparison of lU-pair production with fermion pair production considered in 
the literature, and our conclusions. 

2 Z' Models 

The presence of an additional neutral gauge boson (Z 1 ) is anticipated in many 
extensions of the SM. Some grand unified theories (GUT) like those based on 
E%, which contain the gauge group of SU(3)c x SU(2) L x U(l)y x U(l)' clearly 
have additional 17(1)' symmetries (for a recent review, see, e.g. ref. [13]). These 
could also arise in superstring theories. There are many other extensions of 
the SM with dynamical symmetry breaking, Little Higgs models (LHM), LRSM 
and ALRSM also with extended gauge sectors. In models with large extra 
dimensions, Kaluza-Klein excitations of the SM gauge bosons propagating in 
the bulk manifest as extra gauge bosons in four dimensions. For our purposes 
we confine our attention to certain GUT models based on Eq and LRSM and 
ALRSM. We will study the presence of an extra U(l)' symmetry present in 
addition to the SM gauge symmetries, arising in the candidate models mentioned 
above. In some of these models there could be more than one neutral gauge 
bosons along with possible presence of heavy charged gauge bosons. We assume 
that any such additional gauge bosons decouple from the particle spectrum 
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under study, and therefore will be ignored in this study. 

The new gauge boson Z' could mix with the SM gauge boson to give the 
physical eigenstates. As explained later in this section, the W-pair production 
in e + e~ collisions has the advantage of directly probing the mixing unlike, for 
example, the fermion pair production process. This is because, the W does not 
interact directly with the Z' . While such mixing is highly constrained by the 
current experimental constraints, it may still be possible to probe its effect in 
a high energy, high luminosity machine like the ILC. With very high statistics 
expected for W-pair production at ILC, this process has the potential to probe 
even very small mixing permitted. 

Let us now turn to some general features of the scenario. Here we closely 
follow the discussion in ref. [13l [14] . With one additional U(l)' symmetry, the 
mass term of the neutral gauge bosons may be written as, 



lt ss = Uz^ z">)[ f _ ][ %! ) (i) 




Diagonalizing the above mass matrix the mass term is presented in terms of the 
physical boson fields as 

LT SS = \mI z»z^ + \ul Z$Z 2 ^, (2) 

where we identify the lighter Z\ as the observed Z-boson with M\ — 91.19 GeV, 
and the Z 2 as its heavier counterpart. In terms of the SM gauge boson, Z and 
the new gauge bosons, Z' , we may write the physical states as 



Zi \ _ f cos9 sin (9 \ f Z° 
Z 2 ) ~ { -sm9 cose ) \ Z 1 



(3) 



The mixing angle, 6 is related to the diagonalization of the mass matrix, and 
can be expressed in terms of the physical masses and the SM mass parameter 
as 



ML - M? 

tan^ 6 = - J - * . 4 
Mf-M§ D v ' 



For our phenomenological analysis it is more convenient to re-parametrize the 
above by defining a mass difference AM = Mi — M^o . Rearranging Eq. 0] the 
mass of the heavier gauge boson takes the form 

2 (l+tan 2 0)(M 1 -AM) 2 -M 1 2 

= to&e ' (5) 

Thus, in our study we consider e and AM as the two independent parameters 
of the mixing. The importance of mixing in VF-pair production at ILC is clearly 
visible with the fact that the new gauge sector does not interact with the SM 
gauge sector directly, and therefore the W boson does not couple directly to Z' at 
tree level. Its couplings to the mass eigenstates Z 2 arises through mixing, along 
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with a corresponding weakening in its coupling to the lighter mass eigenstate, 
Z\. These couplings are given by 

gwwzt = 9wwz° cos 0, 9wwz 2 = 9wwz° sm Oi (6) 

where gwwz° is the SM coupling. Experimental constraints on the mixing 
limits the value of 9 to be not larger than a few times 10 -3 . At the same time, 
when the Higgs structure of the model is known, one may compute the mixing 
angle. In such a case, mixing angle can be expressed as 

[5~ Ml 
9 = C J-Xsm9 w j^, (7) 

Zi 

where A is a parameter of order unity, C is function of the VEV's of the Higgs 
fields and their U(l)' charges and 9\y is the usual Weinberg mixing angle. In 
Table [1] we present an illustrative case of E 6 models as discussed in Ref.[5]. The 
table also gives the mixing angle, 9 and AM in each case corresponding to a 
representative value of Mz 2 = 1 TeV. 



Model 


C (range) 


9 


AM (MeV) 






3_ _2_ 






[-.0037, .0025] 


[75, 33] 


Mi>) 




-vW 


2 
5 




[-.0032, .0032] 


56 


E 6 (V) 




' L_ 




[-.0010, .0040] 


[6,89] 


LRSM 




1_ l~z 

Oil! V 5 ' 






[-.0019, .0048] 


[20, 124] 



Table 1: Mixing angle (9) and AM corresponding to a Z2 of mass 1 TeV in 
different models considered. The parameter in the left-right symmetric model 
(LRSM) takes a value, am = \J \ — 2 sin 2 9w / sm9w- 



We now turn our attention to the Zee coupling. The fermion couplings 
are highly model dependent. In the following we will very briefly describe this 
coupling in the models considered. A detailed analysis can be found in the 
literature including |13j . The neutral current interactions of electrons with the 
gauge bosons are given by the Lagrangian term, 

C NC = J? m A^ + J$ M Z° + J'"Z;, (8) 

where J£ m is the electromagnetic current with which the photon interact, Jg M 
is the current with which the SM neutral gauge boson, Z Q interact, and J' M is 
the current with which the new neutral gauge boson, Z' interact. In terms of 
the projection operators Pl.r = (1 T 7 5 )/2 the currents take the form 

Jt = {qIlPl + 9!rPr)A ■ (9) 
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The SM couplings involved in Jg M are expressed in terms of the Weinberg 
mixing angle, Ow as 



9l 



sin t>w cos t)w 



sin 2 0w), 



f e sin t>w 

9r = COS0 W ' 



(10) 



Fermion couplings corresponding to the current, J' M in the cases of models 
considered in our study are tabulated in Table [2j 



Model 



g'l 



E 6 (x) 

E 6 (v) 
LRSM 

ALRSM 



COS 8\v 
TOe 



12 cos 9 w 

e 

6 cos &w 
e 1 

COS 9\y 

— a e - 2n —(-h +sin 2 Ow) 

cos 9w sin Ow ollb. v 2 vv * 



cos % 



24 cos 

12 cos Ow 
e 

3 cos Ow 

( 1 _ OiJ 

V2q L7? . 2 



e . 5 fl — ( — h + § sin' 2 

sin^ W V 2 2 



) 



Table 2: Z'ee couplings (Eq. in different models considered. The parameter 
otLB. = yl — 2 sin 2 6w / sin^w- 



In the next section, we will compare the effects found in the above models, 
with those arising in the Littlest Higgs scenario whose collider signature at the 
ILC with polarized beams was recently considered in ref. [10]. Briefly stated, in 
this scenario, the Higgs fields are considered to be the Nambu-Goldstone Bosons 
(NGB) of the non-linear realization of some global symmetry breaking. 

In our numerical analysis for LHM, there are two free parameters / and Oh- 
As argued by [15] . precision electroweak measurements restrict the parameters 
to be / ~ 1 TcV and 0.1 < cos Oh < 0.9. In our study we consider a value of 
/ = 1 TeV and cos Oh = 0.3 satisfying these restrictions. 

3 Analyses of e + e~ -)• WW 

In this section we present the results of our numerical analysis to probe the 
Z' models through the process e + e~ — > WW at the ILC. This process gets an 
additional s-channel contribution through the exchange of the heavier gauge 
bosons, Z 2 - Only the SM component of Z 2 interacts with W, and therefore this 
additional s-channel contribution is not present in the absence of mixing. At the 
same time, Z\ becomes the SM gauge boson, Z° with the standard couplings, 
when there is no mixing. The t-channel with ^-exchange is not affected by new 
physics here as W interactions with SM fermions proceed through standard 
couplings. Thus the e + e~ — > WW process directly probes the Z° — Z 1 mixing, 
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unlike processes like e + e~ — > //, which gets contribution from the presence of 
Z' even when there is no mixing. 

Thus for the process e + e~ — » W + W~, we have an s-channel process with 
the exchange of the heavy neutral gauge boson, Z 2 , in addition to the standard 
channels as shown in Fig. [T] 



(WWW 

w+ 

WWW 




Figure 1: Feynman diagrams contributing to the process e + e — > W + W in 
the Z' models. 



The three-point gauge couplings involving WW are given by: 
FWr(fe)F(fe) = igvww [g^{ki - k 2 )P + g v P{k 2 - + g^(k 3 - hi)"] , 

where all the momenta are considered outflowing, and V = 7, Z\, Z 2 . 

Direct search results from pp collisions at TeVatron as well as precision 
electroweak analysis constrain the parameters discussed in the previous section 
|5 . In most cases Z 2 mass slightly smaller than 1 TeV is permitted. This is 
translated into allowed range of couplings in a specific version of Eq model as 
presented in Table [U In our numerical studies we consider a conservative value 
of 0.003 and 100 MeV for the mixing angle and the parameter AM, respectively. 
The above parameter choices correspond to a Z 2 of mass ~ 1400 GeV . 

For a non-exhaustive list of phcnomcnological studies of Z' in the context of 
LHC as well as ILC, see reference in [16l [17l HH [19] . Previous phenomenological 
studies of e + e~ — > WW process in the context of Z' have considered different 
observables at the production level [21 [H] , and have obtained the reach of ILC 
in probing the model. In this work our main focus is on various observables con- 
structed with the decay products of the W's produced. One obvious advantage 
is that these observables do not require full reconstruction of the WW events, 
unlike in the former case. In this section, we will first report our analysis of 
the process at the production level with a study of the total cross section, and 
subsequently analyze different decay distributions and other observables. 

3.1 The total cross section 

We compute the total cross section incorporating beam polarization using the 
helicity amplitudes given in ref. [3] with the new couplings and with the added 
contribution due to the exchange of Z 2 . With beam polarization, in general, 
the polarized cross section may be expressed as: 
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a(e+e- ^W+W-) = \ [(1 + P e -).(1 - P e+ )a RL 

+ (l-P e -).(l + P e+ )a LR ], (11) 

where a RL = a(ete R -> W + W~) and cr^ = o"(e^e£ -> W+W - ), with 
e L.fl representing the left- and right-polarized electrons (and positrons), respec- 
tively. The degree of polarization is defined as: P e = (Nr — Nl)/(Nr + Nl), 
where Nl,r denote the number of left-polarized and right-polarized electrons 
(and positrons), respectively. More than 80% of electron beam polarization and 
large positron beam polarization are expected to be achieved at ILC. In our 
analysis we consider both the ideal possibility of 100% polarization along with 
the realistic polarization of the beams that will be achieved at ILC. 

The models considered here are insensitive to the total cross section even 
at higher centre of mass energies. The inclusion of beam polarizations has no 
significant effect. In Table EJ we present the total production cross section for 
SM and Z models considering both unpolarized and polarized beams. It can 
be seen that at -^=500 GeV, the deviation for different Z' models from the 
SM using unpolarized beams is about 0.05 - 0.3%. The deviation is slightly 
increased to about 0.7% at -^5=800 GeV. Switching on the polarization with 
left-handed electrons and right-handed positrons it is seen that with this choice 
of polarization the cross section is about four times more than the unpolarized 
case. Although the deviation in this case follows the same pattern as the unpo- 
larized one, but due to the larger size of the cross section Z — Z' mixing can be 
studied more effectively. However with right-handed electrons and left-handed 
positrons, due to the absence of the dominant t channel the mixing effect is 
more pronounced, even though the cross section is very small compared to the 
other two cases. This particular combination thus leads to an increased signal 
by background ratio. It can be seen that the deviation is about 24% for different 
E 6 models and 50% for the LRSM and ALRSM model at Vi=800 GeV, with 
100% beam polarization. We have only presented the figure for this particu- 
lar case. In Fig. [5] we present the total production cross-section in the case of 
SM and Eq(x) as this has a maximum deviation compared to other Eq mod- 
els. Moreover, it can be seen that the percent deviation of LRSM and ALRSM 
are equal with this specific choice of beam polarization. LHM is also plotted 
for comparison. Since the gauge structure of LHM is different from the other 
Z' models considered here, it behaves differently from them. It is not much 
constrained compared to the other models as can be seen from Fig. [2J Taking 
into account the polarization which will be achieved at ILC, different E§ models 
show about 4% deviation from the SM. The percent deviation for LRSM and 
ALRSM is about 9% compared to 50% obtained in the ideal case. 

3.2 Double Energy Distribution 

In order to exploit further the process at hand, it is profitable to consider the 
decays of one or both the W's. Let us consider e + e~ — > W + W~ with both 
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p- 

1 e 


P+ 


itiocIgI 


a/s —300 GcV 


a/s —500 GeV 


xfs —800 GeV 


,/s —1000 GcV 






SM 


13.598 


7.208 


3.744 


2.692 








13.618 


7.227 


3.767 


2.724 






R c <ih) 


13.602 


7.212 


3.749 


2.699 


n 




^vl) 


13.612 


7 223 

1 . ZZ(J 


3 762 


2 71 7 






LRSM 


1 3 601 


7.212 


3 750 


2 701 

Zj . 1 VJ _L 






ALRSM 


13.571 


7.183 


3.720 


2.666 






LHM 


1 3 679 


7 347 


3 867 


2 973 






SM 


53.973 


28.716 


14.940 


10.746 








54.042 


28 777 


1 4 994 

14:. £7c74 


1 798 








54 005 


28 743 


1 4 963 


1 767 


1 

j. 


1 

_L 




^4 038 


9R 77 A 


1 4 QQ1 


1 7Q^ 






LRSM 


^4 01 ^ 




1 4 981 

14. C/O _1_ 


1 &00 










9R fi40 




lu.UUu 






LHM 


54.214 


29.205 


15.362 


11.796 






OTV f 

SM 


41.023 


21.826 


11 OCC 

11.355 


8.167 








41.076 


21.871 


11.396 


8.206 








41.048 


21.845 


11.372 


8.183 


-0.9 


0.6 


77" / \ 


41.073 


21.869 


11.393 


8.204 






LRSM 


41.053 


21.850 


11.376 


8.187 






ALRSM 


40.964 


21.767 


11.295 


8.102 






LHM 


41.208 


22.198 


11.676 


8.966 






SM 


0.418 


0.114 


0.037 


0.023 






77" / \ 


A A O T 

vA27 


a 1 on 

0.122 


0.045 


0.030 






77" / M 


0.401 


0.103 


0.028 


0.015 


1 


-1 


77> / ^ 


0.411 


0.110 


0.034 


0.020 






LRSM 


0.390 


0.093 


0.018 


0.003 






A T T") OA J 

ALRSM 


0.390 


0.093 


0.018 


0.003 






LHM 


0.500 


0.182 


0.105 


0.096 






SM 


0.857 


0.374 


0.178 


0.125 






E e (x) 


0.865 


0.381 


0.184 


0.131 








0.845 


0.365 


0.171 


0.119 


0.9 


-0.6 


E 6 (v) 


0.853 


0.371 


0.176 


0.123 






LRSM 


0.837 


0.360 


0.167 


0.115 






ALRSM 


0.835 


0.357 


0.162 


0.109 






LHM 


0.922 


0.431 


0.233 


0.191 



Table 3: Total cross section (in pb) for different models with both polarized and 

unpolarized beams. The parameter used for LHM is /=1 TeV and c=0.3, and 

for the Z' models = 0.003 and AM = 100 MeV 

9 




Figure 2: Total cross section for W + W production in an e + e collision for 
SM[Red-Solid], ALRSM, LRSM [Magenta-Dashed], LHM [Green-Dotted] and 
Eq(x) [Blue-Dotted] with polarized beams with P e -=1 and P e +=-l. For U(l)' 
models, 9 = 0.003 and AM = 100 MeV are considered. The parameter values 
of f=l TeV and c=0.3 are used for LHM 



VF's decaying into leptons. The differential cross section in this case can be 
expressed as 

da 



dcosO dcos#* d$*_d cos 6* + d<j>+ 
^—B(W- -> l-u l )B(W+ l+^V&Vl, V\ n (12) 



where O is the scattering angle and /3 = — Arri^ / s) is the velocity of the 
W in the centre of mass frame, y/H being the centre of mass energy. The other 
angles, 6Z. and (j)Z. are the polar and azimuthal angles of the lepton/antilepton 
in the rest-mass frame of W T with boost direction of W~ along the 2-axis, 
respectively. The production and decay tensors, Vvh,, T>\, and "D\, are given 
in [3]. The energy of the l T in the centre of mass frame is related to 9*^ in the 
following way: 

Vs 
4 

This allows us to obtain the double energy differential cross section from Eq. [T^] 
as 



E lT =X- (l±/3cos0;), (13) 



da f f 2 \ da , 

■ a cos B i 



dEi-dE t + J \f3jmw J dcos® d cos 6*_ d<p*_d cos 0+ 

(14) 
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Following the notation of Dicus and Kallianpur [T^] , the energies of the leptons 
Ei± are expressed as dimensionless variables X t ± defined as: 

Xi± varies between and 1. In Table |4] we present the percentage of events in 
bins of X t - and Xi+ obtained from 

1 cPa 
a dXi+dXi- 

at y/s — 800 GeV using unpolarized beams. As expected in the case of SM 
the distribution peaks at maximum values of X t - and X t + . This behaviour 
dominates in other models as well. Since the matrix in Table 2] is symmetric 
under interchange of X[- and X[+ (see Ref. [12] for explicit expression), we 
have shown only the upper half. We have numerically checked that the matrix 
is indeed symmetric. 

A combination of the effects result in about 5% deviation in the case of E§(x) 
and ALRSM models, but with a qualitative difference. In the case of Eq(x) it 
is an enhancement, whereas in the ALRSM case there is a reduction. The effect 
is reduced to smaller than 4% in the case of E^rf) model, while in the case of 
Ee(ip) and LRSM models it is a very negligible contribution of about 1%. The 
most sensitive bin in all the cases is the one with 0.2 < (A;- , X t +) < 0.4, with 
about 0.6% of the total events in the case of SM. Such small deviation of a 
few percent is therefore hard to detect even at a high luminosity machine like 
the ILC. On the other hand, the LHM model, while showing similar qualitative 
behaviour, deviates from the SM value by about 34% in the same bin, leaving 
scope of detection at the ILC. 

In the above analysis we have not made any attempt to optimize our results 
by considering different binning options. We expect the qualitative behaviour 
to remain more or less the same even in the optimal case. This is supported by 
Fig. [3] and Fig. 0] where we plot the double energy distribution for SM and the 
deviation from SM expressed as 



d 2 cr 


d 2 a 




dX l -dX l+ 


, , ClX, — C&C,-L 

!■:■>■!< 1 1 ~ 


SM 


6? a 






dX t -dX l + 


SM 





We can see from Fig.[3]as discussed before, in the case of SM the distribution 
peaks at maximum values of X t - and A ; +. 

The use of right-handed electron beam and left-handed positron beam is 
expected to have a much larger sensitivity to all the models. In the case of 
LHM, while the new gauge boson do not contribute, owing to the fact that this 
does not couple to the right-handed electrons, the changed fermionic couplings 
of the SM gauge boson provides substantial effect [TU]. In this case, the t- 
channel contribution is not present, and therefore both the SM as well as new 
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Table 4: Percentage of events at y / s=800 GeV with unpolarized beams in bins 

of X t - and X t + corresponding to different models: a =SM, b—E 6 (x), c=E e (ip), 

d=E 6 (r}), e=LRSM, /=ALRSM and 5 =LHM. The parameters for U(l)' models, 

6»=0.003 and 5M=100 MeV, while for LHM / = 1 TeV and c = 0.3 are used 
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SM 




Figure 3: Double energy distribution for the SM at y/s — 800 GeV with unpo- 
larized beams 

physics contributions show a symmetric behaviour. The difference in the two 
cases is, as expected, a constant shift, either positive or negative. At y/s = 800 
GeV, about 50% deviation is seen for both LRSM and ALRSM models, while 
Ee(ip) has about 40% deviation. The other two models show slightly reduced 
sensitivity with about 30% and 12% deviations in the case of Eq(x) and E 6 (r]) 
models respectively. 

3.3 Single Energy Distribution 

Energy distribution of the secondary lepton obtained by integrating the double 
energy distribution in Eqn. 1141 over £7+, is another observable that might give 



E 6 ,v , LHM 




Figure 4: Double energy distribution showing the percentage deviation for the 
two models (a) E 6 (x) and (b) LHM from the SM at yfs = 800 GeV with unpo- 
larized beams. In the case of E 6 (x) = 0.003 and AM = 100 MeV, and in the 
case of LHM / = 1 TeV and c = 0.3 are used 
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a handle on the new effects. In Figs. |5] we plot the energy distribution of 
the lepton for the case of y/s = 800 GeV for both ideal and realistic degrees 
of beam polarizations. The signs of the beam polarization are chosen so as to 
essentially switch off the SM t-channel contribution, so as to enhance the effects 
of the new physics. There is no appreciable deviation in the case of unpolarized 
beams (and in the case of left-handed electron beams). We notice that the 
LRSM and ALRSM cases behave qualitatively differently compared to the E§ 
models and the LHM. While in the former case there is a reducing effect for the 
entire range of Xi-, the latter has an increasing effect. This could be used as 
a discriminating factor between the left-right symmetric models and the other 
models used. 



0.(115 




Figure 5: Energy distribution of the sec- 
ondary leptons. with polarized beams with 
P e -=1 and P e +=-l at y/s = 800 GeV. 
The parameters used are AM=100 MeV 
and 0=0.003 for U(l)' type of models, and 
/ = 1 TeV and c = 0.3 in the case of 
LHM. Different curves correspond to SM 
(Red-Solid), LHM (Green-Dotted), E 6 {x) 
(Blue-Dotted), and ALRSM and LRSM 
(Magenta-Dashed) . 




0.0 0.2 0.4 0.6 0.8 1.0 

x r 

Figure 6: Energy distribution of the sec- 
ondary leptons. with polarized beams with 
P e -=0.9 and P e+ =-0.6 at y/s = 800 GeV. 
The parameters used are AM=100 MeV 
and 0=0.003 for U(l)' type of models, and 
/ = 1 TeV and c = 0.3 in the case of 
LHM. Different curves correspond to SM 
(Red-Solid), LHM (Green-Dotted), E e {x) 
(Blue-Dotted), and ALRSM and LRSM 
(Magenta-Dashed) . 



Note that in the case of realistic degrees of polarization since the neutrino 
t- channel effects are present, the panels comparing ideal and realistic degrees 
of freedom appear quite different as functions of AT;-. 

3.4 Angular Spectrum of the Secondary Lepton 

We next consider the angular distribution of one of the secondary leptons. The 
way of calculating the angular distribution is done in our earlier work |10j . 
which we follow here. The angular distributions for different polarization com- 
binations in case of different models is calculated for yfs =800 GeV. As in the 
earlier cases, the case of unpolarized beams and the case of left-handed electron 
beams are not significant in the case of angular distribution as well. But the 
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case of right-handed polarization of electron beam along with the left handed 
polarization of the positron beam provides much better discrimination. This 
distribution is shown in the Figs. [7J [H] for both the ideal and realistic degrees 
of beam polarizations. Qualitatively, the picture as regards discrimination be- 
tween the models, remains more or less the same as in the case of single energy 
distribution, except that there is a small difference between the two cases of 
left-right symmetric models considered. In the best case, the ALRSM model 
shows a deviation of around 45% for the parameter set used. 




Figure 7: Angular distribution of one 
of the secondary leptons at ^=800 
GeV using polarized beams with P e -=1 
and P e +~ -1 in the case of SM (Red- 
Solid), LHM (Green-Dotted), E 6 (x) (Blue- 
Dotted), ALRSM (Yellow-DotDashed) and 
LRSM (Magenta-Dashed). The parameters 
used are AM=100 MeV and 9 = 0.003 in 
the case of Eg and LR models, and / = 1 
TeV and c = 0.3 in the case of LHM. 



Figure 8: Angular distribution of one 
of the secondary leptons at ^=800 GeV 
using polarized beams with P e _ =0.9 and 
P e += -0.6 in the case of SM (Red- 
Solid), LHM (Green-Dotted), E 6 (x) (Blue- 
Dotted), ALRSM (Yellow-DotDashed) and 
LRSM (Magenta-Dashed) . The parameters 
used are AM=100 MeV and 9 = 0.003 in 
the case of E^ and LR models, and / — 1 
TeV and c = 0.3 in the case of LHM.. 



Fig. [5] shows the forward-backward asymmetric behaviour arising from the 
Zee coupling, which is different for different models considered. With this obser- 
vation, we may obtain the fraction of leptons emitted in the backward direction, 
which may be defined as 

J-i (da/dcosOi) dcos9i 

fback — ~T[ • 

J-i (da/dcos9i) dcos9i 

This fraction is a useful quantity to consider in the case of unpolarized and 
left-polarized electron beams as well. In Table [5] we present these fractions for 
y/s = 500 GeV, 800 GeV and 1000 GeV. The result shows with unpolarized 
beams about 3-4 % of the leptons are emitted in the backward hemisphere. The 
deviation is about 3% in going from SM to Eq(x) and ALRSM which are more 
sensitive compared to other models. At ^/s=800 GeV, the deviation becomes 
more significant and is about 8%. However LHM being more sensitive shows 
~ 36 % deviation at y / s=500 GeV which further increases with energy. This 
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above results are slightly increased by switching on the beam polarization to 
left-handed electrons and right-handed positrons. The departure from ideal 
degrees of polarization does not play a significant role at the level of significant 
places retained in the tables for this configuration. 

We see from Table [S] that for the case of polarization with about 90% of 
right handed electrons and 60% of left handed positrons, the deviation increases. 
Owing to the near absence of the dominant t- channel neutrino diagram with this 
polarization configuration, the new physics contribution shows up in the fraction 
of the leptons emitted in the backward direction. Thus at \/s= 800 GeV for 
Ee(tf>) and E 6 (x) there is about 15% deviation from SM, whereas E 6 (r]) is more 
constrained with only 5% deviation. The symmetric models ALRSM and LRSM 
are more sensitive with about 30% deviation. 

In the case of completely right handed electrons and left handed positrons, 
from the same table it may be seen that all the models give rise to the same 
prediction as the SM. This is a consequence of the near symmetric behaviour 
in the angular distributions as shown in Fig. [7] Due to this property the new 
physics effects in fback and Afb are essentially wiped out. It is important 
to note that this feature persists with completely right-polarized electron and 
partially left-polarized positrons, and vice versa. 

Another useful observable related to the angular asymmetry is the forward- 
backward asymmetry defined as 



The values for this asymmetry is also presented in Table [5] for different models at 
three different collider energies. This asymmetry shows a similar behaviour as 
that of the fraction of the leptons emitted in the backward direction as regards 
the beam polarization configurations and discrimination between the models. 
In fact, the two observables are intimately related: Afb = 2fback — 1- How- 
ever, since they are experimentally realized differently, it may be appropriate to 
consider both of them. 

3.5 Left-Right Asymmetry 

The left-right asymmetry is another important observable to be considered for 
studying the new models. We define the left-right asymmetry in the differential 
cross section as: 



where 9 is the W scattering angle. 

The Z' of the U(l)' models considered here couples to both left- and right- 
handed fermions, but with varying relative couplings. Thus, one would expect 
appreciable change in the asymmetry between the left- and right-polarized cross 



Afb 




(16) 



J_ 1 (da jd cos 9 '/) dcos9i 



di f f _ da (e^e L )/d cos 6 - da (e^ e R ) / d cos 9 
da (e^e£)/d cos 9 + da (e£e^)/d cos 9' 



(17) 
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y/s =500GcV 


y/s =800GcV 


y/s =1000GeV 


Pe- 


Pe+ 


Model 


f, 

J back 


afb 


f, 

J back 


afb 


f, 

J back 


afb 






SM 


0.0346 


-0.9308 


0.0238 


-0.9524 


0.0208 


-0.9585 






^6 (A) 


0.0356 


-0.9288 


0.0257 


-0.9487 


0.0234 


-0.9532 








0.0347 


-0.9305 


0.0242 


-0.9517 


0.0213 


-0.9574 


o 


o 




0.0353 


-0.9293 


0.0252 


-0.9496 


0.0228 


-0.9545 






LRSM 


0.0347 


-0.9306 


0.0241 


-0.9517 


0.0213 


-0.9573 






ALRSM 


0.0330 


-0.9339 


0.0208 


-0.9585 


0.0161 


-0.9677 






LHM 


0.0474 


-0.9053 


0.0474 


-0.9052 


0.0212 


-0.9574 






SM 


0.0324 


-0.9353 


0.0225 


-0.9550 


0.0197 


-0.9606 






E«(V) 
J^e \a) 


t~\ noon 

0.0332 


-0.9335 


0.0241 


-0.9518 


0.0219 


-0.9561 








0.0327 


-0.9345 


0.0232 


-0.9537 


0.0206 


-0.9588 


-0.9 


0.6 


E 6 (v) 


0.0332 


-0.9336 


0.0240 


-0.9520 


0.0218 


-0.9563 






LRSM 


0.0328 


-0.9343 


0.0233 


-0.9533 


0.0209 


-0.9583 






ALRSM 


0.0312 


-0.9376 


0.0201 


-0.9598 


0.0160 


-0.9681 






LHM 


0.0440 


-0.9120 


0.0441 


-0.9118 


0.0167 


-0.9666 






SM 


0.1618 


-0.6765 


0.1061 


-0.7878 


0.0909 


-0.8183 






E 6 ( X ) 


0.1694 


-0.6612 


0.1208 


-0.7583 


0.1116 


-0.7788 






E e W 


0.1518 


-0.6964 


0.0889 


-0.8222 


0.0680 


-0.8639 


0.9 


-0.6 


Eo(v) 


0.1587 


-0.6826 


0.1013 


-0.7975 


0.0845 


-0.8310 






LRSM 


0.1457 


-0.7087 


0.0786 


-0.8427 


0.0552 


-0.8896 






ALRSM 


0.1426 


-0.7148 


0.0648 


-0.8703 


0.0268 


-0.9463 






LHM 


0.2174 


-0.5653 


0.2120 


-0.5760 


0.2251 


-0.5498 


1 


-1 


All models 


0.5834 


0.1668 


0.5392 


0.0785 


0.5264 


0.0528 



Table 5: Fraction of leptons emitted in the backward direction, and the forward- 
backward asymmetry for all models for unpolarized and polarized beams with 
a parameter choice of 0=0.003 and AM= 0.1 GeV. The parameters for LHM 
are /=1 TeV and c=0.3 
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sections. But compared to these the new gauge boson in the LHM has the 
peculiar property of coupling only to left-handed fermions as mentioned earlier. 

Fig. Ufa) shows the LR asymmetry for -^=800 GeV. Even at low energies the 
deviation becomes apparent. Notice also that the LHM now pair with the left- 
right symmetric models unlike in the case of energy and angular distributions 
shown in Figures [5] and [JJ providing us another tool for discriminating different 
models from each other. 




200 400 600 800 1000 



cos9 VT(GeV) 



Figure 9: (a) Differential Left-right asymmetry as a function of the scattering 
angle for different models at y/s — 800 GeV, and (b) the integrated left-right 
asymmetry as a function of the centre of mass energy for different models: 
SM (Red-Solid), LHM(Green-Dotted), E 6 (x) (Blue-Dotted), LRSM(Magenta- 
Dashed) and ALRSM (Ycllow-DotDashed) . Parameters used are 9 = 0.003 and 
AM = 100 MeV for the E 6 and LR models, and / = 1 TcV and c = 0.3 for the 
LHM. 



We may go one step further by considering an integral version of this asym- 
metry as better efficiency may be obtained this way, by integrating each of the 
differential cross sections from an opening angle 9o up to an angle it — 9$, for var- 
ious realistic values of 8a to which the data can be integrated without difficulty. 
We define the integrated left-right asymmetry as: 

A = -> ww) - <jjM^R -> ww) 

LR <7e (e+e- L ^WW) + ag (e+e R ^WW) 

where ag stands for JJ^ 6 ° (da/d9) d6. 

This asymmetry, for different parameter models is plotted against the cen- 
tre of mass energy in Fig. |9j&) with 9$ =0. Dominance of the t-channel in the 
VF-pair production establishes a highly forward peaked cross section, whereas 
as is seen from Fig. Efa) the deviations in LR asymmetry grows with the scat- 
tering angle, except for a region near 9 = 180 degrees. Thus, the cut-off angle 
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can be effectively used to increase the deviation in the asymmetry. The results 
are presented in Table [6] for different cut off angles in different models at dif- 
ferent centre of mass energies. An asymmetric cut off, meaning and integrated 
asymmetry between scattering angles, 6*io and it — 6*20, with #20 7^ #10, may 
improve the situation quantitatively. Again, as in the case of other observables, 
our gaol in this report is to demonstrate the viability of identifying deviations 
from the SM case, and further providing tools to disentangle different possible 
models. With this illustrative purpose in mind, we do not attempt to optimize 
the investigation. 

4 Discussion and Conclusions 

In the present work we have considered a class of additional Z' models which are 
of interest to the linear collider community, both at the ILC as well as at CLIC. 
While the masses of these bosons are already required to be significantly high, 
their imprint through mixing with the standard model Z boson is the subject of 
this investigation. By considering popular Eq and left-right symmetric models 
like LRSM and ALRSM we have demonstrated that the new physics signatures 
due to these models can be imprinted only at higher center of mass energies. 
In the LHM model considered in our earlier work [10], Z' had the property of 
coupling only to left-handed fermions. This is in contrast to the other models 
considered here where Z' behaves like the SM Z. Thus compared to LHM, these 
models do not show appreciable deviation at lower center of mass energies. 

While our focus in W-pair production at ILC is on the unambiguous signal 
that it provides for Zsm — mixing, we notice the interesting possibility of 
model discrimination here. For example, let us compare the deviations of some 
of the observables from their SM values. In the case of energy and angular 
momentum distributions notice the qualitatively different behaviour of left-right 
symmetric models compared to E§ models and LHM. Similarly, in the case of 
integrated and differential left-right asymmetries the E$ model has a different 
qualitative behaviour compared to the LHM and left- right symmetric models. 

Note that WWZ. L (i — 1,2) coupling is insensitive to differences between 
Z' models, and therefore it is the coupling of Z' with the initial electron and 
positron that enables any possible model discrimination in the present case. 
Thus, one would imagine that fermion pair production process is better suited 
to distinguish different models. In a recent study, the potential of the ILC to 
discriminate between Z' models through fermion pair production is studied in 
Ref. [TB] • Here the sensitivity is studied by considering cross sections which are 
shown to be sensitive to new physics due to the availability of high statistics at 
the ILC. It is also shown that beam polarization does enhance the sensitivity in 
an essential manner. Similar studies with Drell-Yan dilepton production [17] also 
probe distinguishability of different Z' models at LHC. At the same time, our 
study indicates that VF-pair production at ILC is capable of supplementing the 
fermion pair production process in model discrimination in the case of Z' models. 
More extensive analysis involving a parameter scan, also incorporating a realistic 
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A L r 


Cut off angle 8 


Model 


V^=500GeV 


sjs= 800 GeV 


V«= 1000 GeV 




SM 


0.9921 


0.9950 


0.9958 




E 6 (x) 


0.9914 


0.9936 


0.9930 




E 6 (ip) 


0.9930 


0.9967 


0.9984 





E e (v) 


0.9924 


0.9956 


0.9968 




LRSM 


0.9935 


0.9976 


0.9995 




ALRSM 


0.9935 


0.9976 


0.9995 




LHM 


0.9945 


0.9972 


0.9983 




SM 


0.9853 


0.9872 


0.9877 




E 6 (x) 


0.9842 


0.9836 


0.9800 






0.9870 


0.9916 


0.9955 


15 


E 6 (v) 


0.9860 


0.9889 


0.9909 




LRSM 


0.988 


0.9939 


0.9985 




ALRSM 


0.9879 


0.9938 


0.9985 




LHM 


0.9885 


0.9933 


0.9956 




SM 


0.9742 


0.9765 


0.9771 




E 6 (x) 


0.9723 


0.9703 


0.9638 




E 6 {ip) 


0.9772 


0.9846 


0.9917 


30 


E 6 (v) 


0.9755 


0.9797 


0.9834 




LRSM 


0.979 


0.9889 


0.9973 




ALRSM 


0.9787 


0.9885 


0.9971 




LHM 


0.9808 


0.9890 


0.9928 




SM 


0.9601 


0.9622 


0.9628 




E e (x) 


0.9574 


0.9531 


0.9434 






0.9648 


0.9753 


0.9866 


45 


Ee{v) 


0.9622 


0.9679 


0.9738 




LRSM 


0.9676 


0.9822 


0.9956 




ALRSM 


0.9669 


0.9811 


0.9951 




LHM 


0.9721 


0.9845 


0.9899 




SM 


0.9453 


0.9466 


0.9470 




Ee(x) 


0.9418 


0.9348 


0.9224 




Mi>) 


0.9517 


0.9652 


0.9811 


60 


Ee(v) 


0.9483 


0.9552 


0.9638 




LRSM 


0.9555 


0.9750 


0.9939 




ALRSM 


0.9542 


0.9727 


0.9927 




LHM 


0.9640 


0.9807 


0.9875 



Table 6: Left- right asymmetry for different cut-off angles at selected y/s values 
for different models considered with a parameter choice of 0=0.003 and AM— 
0.1 GeV. The parameters for LHM are f=l TcV and c=0.3 

20 



collider-detector simulation is required to draw conclusions regarding ability of 
VF-pair production in model discrimination at ILC, adapting the procedure in 
Ref. 16 . Alternatively one may contemplate adapting different asymmetries 
of the type we have proposed for fermion pair production as well. Perhaps a 
joint analysis involving both fermion pair production and W^-pair production is 
needed to bring out the full potential of ILC in this regard. 

Our conclusions for our models is as follows: we have studied in detail each 
of the Eg models denoted by \-> i> an d f] and the LRSM and ALRSM models. 
We find that at a centre of mass energy of 800 GeV, all these models show a 
sharp deviation from the SM predictions for energy-energy, single-energy and 
angular correlations of the decay lepton(s) generally for right-handed electron 
and left-handed positron with realistic degrees of polarization. Curiously for 
ideal polarization the decay lepton fraction and Afb are insensitive to new 
physics for the same configuration above, as the t- channel contribution is totally 
absent in this case. The Xi LRSM and ALRSM models are more sensitive than 
the other models even with unpolarized beams. The reason for the absence of 
sensitivity at lower energies is due to the stringent bounds already present on 
the parameters of the models. The FB and LR asymmetries have also been 
studied, with the latter being more sensitive to new physics in all the models. 
These models remain less sensitive than the LHM model studied earlier. Any 
indication of U{1)' type of Z' in VF-pair production at ILC is a clear signal of 
Zsm — Z' mixing. Apart from this, the present study points to the potential 
of this process to supplement the fermion pair production process in model 
discrimination, through a suitable combination of observables considered. 

Overall a strong polarization program at the ILC will significantly enhance 
the diagnostic capability towards additional Z' bosons particularly at higher 
design energies. 
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